Electronic structures of antiperovskite superconductors: MgXNi 3 (X=B,C,N) 
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We have investigated electronic structures of a newly dis- 
covered antiperovskite superconductor MgCNi3 and related 
compounds MgBNi 3 and MgNNi 3 . In MgCNi 3 , a peak of 
very narrow and high density of states is located just below 
Ep, which corresponds to the n* antibonding state of Ni-3d 
and C-2p but with the predominant Ni-3d character. The 
prominent nesting feature is observed in the F-centered elec- 
tron Fermi surface of an octahedron-cage- like shape that orig- 
inates from the 19th band. The estimated superconducting 
parameters based on the simple rigid-ion approximation are 
in reasonable agreement with experiment, suggesting that the 
superconductivity in MgCNi3 is described well by the conven- 
tional phonon mechanism. 
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Recently, He et have discovered a new intermetal- 
lic superconductor MgCNi3 with the transition temper- 
ature T<7 of 8K, which has the antiperovskite structure. 
Because it has a large proportion of Ni per unit cell, 
it is expected that the magnetic fluctuation would be 
important in determining the superconducting behavior. 
This system is reminiscent of another Ni-based suuercon- 
ductors LnNi 2 B 2 C (Ln = Y, Tm, Er, Ho, Lujm The 
band calculations indicate that very large and narrow 
energy peak in the density of states (DOS) is located 
just below the Kwnii energy Ep which has mainly the 
Ni 3d characterBQ. The behavior of the upper critical 
field H C 2 (E) can be well fitted with the conventional BCS 
expressionQ, whereas the zero-bias conductance peak ob- 
served below Tc suggests that MgCNi 3 is likely to be a 
non-s wave superconductortl 

With Cu doping (electron doping) on Ni-site, Tc de- 
creases systematically, but with Co doping (hole doping) , 
the superconductivity disappears abruptly for doping of 
only In the case of Co doping, there is no evi- 

dence that the quenching of superconductivity is related 
to magnetism. Furthermore, it is observed that the su- 
perconductivity of MgC x Ni3 is sensitive to the content of 
C; it disappears between x=0.96 and x=0.90tl 

To understand the mechanism of the superconductiv- 
ity in MgCNi3, we have investigated systematically the 
electronic structures of MgXNi3 (X=B, C and N), which 
have different number of valence electrons, but have the 
similar band structures. Using the linearized muffin-tin 
orbital (LMTO) band method in the local density ap- 
proximation (LDA), we have obtained band structures, 
DOSs, and Fermi surfaces, and discussed the bonding 
characters. Muffin-tin orbitals up to d-states for Mg, C, 
and up to /-states for Ni are included in the LMTO band 
calculations. We have also estimated superconducting 



parameters based on the rigid-ion approximation. We 
have considered cubic MgCNi3 with the lattice constant 
of 3.81 AEl and employed the atomic radii of 3.20, 1.54 
and 2.49 A for Mg, C and Ni, respectively. The same 
structural parameters are used for all MgXNi3. 

MgCNi3 has the cubic antiperovskite structure: Mg at 
(0 0), C (0.5 0.5 0.5), and Ni at (0.5 0.5 0), (0.5 0.5) 
and (0 0.5 0.5). It is called as an antiperovskite structure 
because the transition metals are located at the corners 
of the octahedron, cage in contrast to the ordinary per- 
ovskite structurea. Without C located at the center of 
cubic cell, MgNi3 is a simple ordered intermetallic com- 
pound with fee structure (Cu3^4u-type). Without C, the 
Ni-3e? band of MgNi3 is very narrow, leading to a mag- 
netic ground state with Ni magnetic moment of 0.43/1 
By inserting C, two carbons become nearest neighbors of 
Ni, and thus it is expected that Ni-3d and C-2p electrons 
are strongly hybridized. 

Total and projected local DOS of MgCNi3 are provided 
in Fig. 0. The overall shape of the total DOS is similar to 
those of Hayward et ala and Dugdale and JarlborgQ. The 
peaks near — 7eV and 4eV correspond to a bonding and 
antibonding states, respectively, of Ni-3c? and C-2p states. 
On the other hand, the peaks near — 4eV correspond to 
ir bonding states of Ni-3d and C-2p. The tt* antibonding 
states are located just below Ep, yielding the high DOS 
at E F (N(E F )) of 5.34 [states/eV]. The contribution of 
Ni-3d states to the DOS at Ep is as much as 76%. Small 
amount of C-2p states are hybridized with Ni-3d states. 
Because the peak just below Ep is very high and narrow, 
this system is expected to be unstable by small perturba- 
tion. The peak is located ~ 60meV below Ep, and about 
0.5 electrons are occupied between the peak and Ep. The 
stoner parameter S, defined as S = N(Ep)Ixc with Ixc 
denoting the intra-atomic exchange-correlation integral, 
is 0.64. Indeed the filling of holes in a rigid band scheme 
produces the magnetic instability. That is, replacing Ni 
by a virtual atom with atomic number 27.93 (correspond- 
ing to CO|-2% doping) yields the stoner parameter larger 
than l.odia. 

The band structure of MgCNi3 along the symmetry 
line of the simple cubic Brillouin zone is shown in Fig. |^. 
The band near — 12eV corresponds to C-2s states, while 
the dispersive bands in the range of — 7eV to — 4eV are 
due to G-2p states. Only two bands (the 18th and 19th 
bands), which have mainly Ni-3c? character, cut the Fermi 
level. These two bands are confined between — 0.5eV and 
l.OeV, and available states in these bands are about four. 
Hence, in the rigid band scheme, small electron or hole 
doping will produce the carriers with the 18th and the 
19th band character. The 18th band is relatively flatter 
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than the 19th band, and so the Ni-3d character is stronger 
in the 18th band. 

Figure [3] presents the Fermi surfaces of the 18th (a) 
and the 19th bands (b) in the simple cubic Brillouin zone. 
The 18th band gives rise to a clover-like hole surface cen- 
tered at X of each cubic face and small hole pockets along 
the (111) directions. With increasing the Fermi level 
(electron doping), the areas of these hole Fermi surfaces 
decrease, reducing N(Ep). In opposite, with decreasing 
the Fermi level (hole doping), the Fermi surface areas in- 
crease, enhancing N(Ep). When decreasing the Fermi 
level further below the DOS peak, the crossing of each 
hole surface occurs, converting it to electron surface and 
decreasing Fermi surface. On the other hand, the 19th 
band yields an octahedron-cage-like electron surface cen- 
tered at r and additional narrow electron surfaces along 
the Brillouin zone edges (Fig. ||(b)). As compared to 
the case of the 18th band, the 19th band shows a rather 
small change in the Fermi surface topology with varying 
the Fermi level position due to its more dispersive band 
character. The Fermi surface topology obtained in the 
present study-iis qualitatively similar to those of Dugdale 
and Jarlborgp. 

The most notable in Fig. |^(b) is the prominent nesting 
feature along the (110) direction observed in the ab plane 
of the T-centered Fermi surface with the octahedron- 
cage-like shape. In fact, this is contrary to the report by 
Dugdale and Jarlborgp who have not observed the obvi- 
ous nesting feature in the Fermi surface of the 19th band. 
As mentioned above, the detailed shapes of the Fermi 
surfaces are very sensitive to the position of the Fermi 
level, because of the very sharp DOS peak near Ep. Pre- 
sumably the difference between two results arises from 
the different band parameters employed in the LMTO 
band calculations, such as atomic radii, number of k- 
points, energy parameters, and so on. In any case, the 
present result reveals that the system is in the vicinity 
of the Fermi surface nesting, if not complete in undoped 
MgCNi 3 . It is well known that the system with the Fermi 
surface nesting can be strongly correlated with various in- 
stabilities: the structural transition (the charge density 
wave instability) or the spin density wave instability. It 
is thus expected that MgCNi 3 may be susceptible to one 
of the above instabilities. However, until now, any ev- 
idence otmagnetic or structure transition has not been 
reportecoEi This aspect remains to be resolved. 

We have seen that the DOS peak near Ep is pro- 
duced by the hybridization of Ni-3d and C-2p states. In 
MgCNi 3 , the bands from — 0.5eV to l.OeV are almost 
half-filled: two electron states out of four available states 
are occupied. To explore the doping effect with varying 
the number of valence electrons, we have investigated 
electronic structures of MgXNi 3 (X= B, C, N) . Figure 
|] provides the DOSs for MgXNi 3 . It is seen that the 
effect of changing X is mainly a variance of the Fermi 
level with respect to the DOS peak. The shape of DOS 
is perturbed a little, which indicates that the rigid band 
scheme would work well in this system. 



The B-2p state in MgBNi 3 is located higher in energy 
than the C-2p state of MgCNi 3 , and so the hybridization 
with Ni-3d is stronger. Hence the band is more dispersive 
and accordingly the DOS peak becomes smeared. Al- 
though the Fermi level in MgBNi 3 is located very close 
to the DOS peak, the DOS at E F , 4.79 [states/eV], is 
comparable to that of MgCNi 3 (5.34 [states/eV]) (Table 
||). Hence the magnetic instability does not occur either 
in MgBNi 3 . As described before, the effective hole dop- 
ing in MgBNi 3 converts a hole Fermi surface of the 18th 
band to an electron Fermi surface, and an electron sur- 
face of the 19th band is reduced smaller. In MgNNi 3) the 
N-2p state is located a bit lower in energy than the C-2p 
state of MgCNi 3 , yielding reduced bandwidths of both 
the 18th and the 19th band. By the effective electron 
doping in MgNNi 3 the DOS at E F is reduced to 3.63 
[states/eV] (Table @). The contribution of the 18th band 
to the DOS at Ep is almost negligible and the Fermi sur- 
face of the 19th band, which gives the main contribution 
to the DOS at Ep, is changed to a hole surface. 

We have explored superconducting properties of 
MgXNi 3 based on the rigid-ion approximation^. We 
have estimated the superconducting parameter rj a = 
N(Ep)(I 2 ), where (I 2 ) is the average electron-ion inter- 
action matrix element for the a-th ion. Table || provides 
the calculated r\ a for each MgXNi 3 . It is seen that the 
contribution of Ni-3d states to the superconductivity is 
most important and r\ Ni is the largest for MgCNi 3 . This 
is consistent with the observed trend that both the elec- 
tron and hole dopings on MgCNi 3 suppress the supercon- 
ductivity. By increasing the atomic number from B, C to 
N, the contribution of X-2p states increases, while that 
of Ni-3d increases first and then decreases. Due to light 
ionic masses of X, even the small increase in rjx affects 
the superconducting property substantially. Therefore, 
effectively electron doped system MgNNi 3 , once synthe- 
sized successfully in the antiperovskite structure, would 
have comparable Tc to MgCNi 3 . 

One can evaluate the electron-phonon coupling con- 
stant X p h by using the McMillan's formula X p h = 
J2 a Va/Maluj^), where M a is aa, ionic mass and (uj^) 
is the relevant phonon frequencytJ. Since there has been 
no information on the relevant phonons, we instead use 
the average phonon frequency (ui 2 ) ~ <d 2 D /2, where 0^> is 
the Debye temperature. However, even the value of ©£> 
is not available. Albeit vera crude, one can estimate ©d 
from the specific heat dataEJ: ~ 300 K from the graph of 
C/T vs. T 2 . Using these informations for MgCNi 3 , one 
obtains X p h — 1.56, and then the McMillan's Tc formula 
with an effective electron-electron interaction parameter 
/i* = 0.13 gives rise to Tc = 23 K. These values seem 
to be too large, as compared to experimental Tc and 
the estimated X p h ~ 0.8 from the specific heat dataEJ. 
Note, however, that X p t strongly depends on the choice 
of the Debye temperature. As shown in Table [n], with 
a choice of larger ©# = 400 K, one obtains X p h = 0.77 
and Tc = 11K, which are in reasonable agreement with 
experiment. This suggests that the superconductivity in 
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this system can be described by the conventional phonon 
mechanism. For more precise estimations of X p h and Tq, 
detailed information on the phonon spectra is prerequi- 
site. 

In conclusion, we have investigated electronic struc- 
tures of the non-oxide antiperovskite superconductor 
MgCNi3. The it* antibonding state of Ni-3d and C-2p is 
formed near Ep with the mainly Ni-3d character. Fermi 
surfaces are composed of two bands. The topology of hole 
Fermi surfaces coming from the 18th band is sensitively 
modified by the variance of the Fermi level position. The 
electron surface of the 19th band with an octahedron- 
cage shape tends to induce the Fermi surface nesting. By 
comparison of DOSs for MgXNi 3 (X=B,C,N), the dop- 
ing effects are discussed. The estimation of X p h and T c 
based on the rigid-ion approximation suggests that the 
superconductivity of MgCNi3 is described well by the 
conventional phonon mechanism. 

Acknowledgments — This work was supported by the 
KOSEF through the eSSC at POSTECH and in part by 
the BK21 Project. Helpful discussions with N.H. Hur are 
greatly appreciated. 



1 T. He, Q. Huang, A. P. Ramirez, Y. Wang, K. A. Reran, 
N. Rogado, M. A. Hayward, M. K. Haas, J. S. Slusky, K. 
Inumara, H. W. Zandbergen, N. P. Ong and R. J. Cava, 
Nature (London) 411, 54 (2001). 

2 R. Nagarajan, Chandan Mazumdar, Zakir Hossain, S. K. 
Dhar, K. V. Gopalakrishnan, L. C. Gupta, C. Godart, B. 
D. Padalia, and R. Vijayaraghavan, Phys. Rev. Lett. 72, 
274 (1994). 

3 R. J. Cava, H. Takagi, H. W. Zandbergen, J. J. Krajewski, 
W. F. Peck Jr, T. Siegrist, B. Batlogg, R. B. can Dover, 
R. J. Felder, K. Mizuhashi, J. O. Lee, H. Eisaki and S. 
Uchida, Nature (London) 367, 252 (1994). 

4 W.E. Pickett and D.J. Singh, Phys. Rev. Lett. 72, 3702 
(1994): L. F. Mattheiss, Phys. Rev. B 49, 13279 (1994); J. 
I. Lee, T. S. Zhao, I. G. Kim, B. I. Min, and S. J. Youn, 
Phys. Rev. B 50, 4030 (1994). 

5 M. A. Hayward, M. K. Haas, A. P. Ramirez, T. He, K. 



A. Regan, N. Rogado, K. Inumara and R. J. Cava, cond 



mat/0104541, (2001) 



b S.B. Dugdale and T. Jarlborg, |cond-mat/010534S| , (2001). 

7 S. Y. Li, R. Fan, X. H. Chen, C. H. Wang, W. Q. Mo, K. 
Q. Ruan, Y. M. Xi ong, X. G. Luo, H T. Zhang, L. Li, Z. 
Sun and L. Z. Cao, |cond-mat/0104554 (2001). 

8 Z. Q. Mao, M. M. Rosario, K. Nelson, K. Wu, I. G. Deac, 
P. Schiffer and Y. Liu, T. He, K. A. Regan and R. J. Cava, 



cond-mat/0105280, (2001) 



D.A. Papaconstantopoulos and W.E. Pickett, Phys. Rev. B 
45, 4008 (1992); P.R. Vansant, PE. Van Camp, V.E. Van 
Doren, and J.L. Martins, Phys. Rev. B 57, 7615 (1998); 
A.L. Ivanovskii, R.F. Sabiryanov, and A.N. Skazkin, Fiz. 
Tverd. Tela, 40, 1667 (1998) (Phys. Solid State 40 1516, 



(1998)); W.S. Kim, E.O. Chi, J.C. Kim, H.S. Choi, and 
N.H. Hur (private communications). 

10 S. H. Shim and B. I. Min (unpublished). 

11 Q. Huang, T. He, K. A. Regan, N. Ro gado, M. Hayward 
M. K. Haas, K. Inumara and R. J. Cava, cond-mat/0105240 
(2001). 

12 G. D. Gaspari and B. L. Gyorffy, Phys. Rev. Lett. 28, 801 
(1972). 

13 W. L. McMillan, Physical Review, 167, 331 (1967). 




-14 -12 -10 -8 -6 -4 -2 2 4 6 
ENERGY (eV) 

FIG. 1. Total and projected local DOS of MgCNi 3 . 
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FIG. 3. Fermi surfaces of MgCNi 3 in the simple cubic 
Brillouin zone coming from the 18th band (a) and the 19th 
band (b) 



FIG. 2. Band structure of MgCNi3 along the symmetry 
lines of the simple cubic Brillouin zone 
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FIG. 4. Density of states of MgXNi 3 (X = B, C, N) 



TABLE I. Total and partial DOSs at E F (in states/eV) for 
MgXNi 3 (X = B,C,N). 



Ni 



Mg 



iVx 



MgBNi 3 0.38 0.18 1.41 4.79 

MgCNi 3 0.22 0.42 1.57 5.34 
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MgNNi 3 0.11 0.34 1.06 3.63 



TABLE II. Comparison of rj (in eV/A 2 ) and \ ph for 
e D =300K and 400K. 







??x 


??Ni 


\ ph (300K) 


\ ph (400K) 


MgBNi 3 


0.00 


0.22 


0.67 


0.67 


0.38 


MgCNi 3 


0.00 


0.48 


1.36 


1.36 


0.77 


MgNNi 3 


0.00 


0.58 


0.87 


1.07 


0.60 



